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Abstract. Deep CCD photometry in the uvby— Hp intermediate-band system is presented for the cluster NGC 2548 
(M 48). A complete membership analysis based on astrometric and photometric criteria is applied. The photo- 
metric analysis of a selected sample of stars yields a reddening value of E(b — y) = 0.06±0.03, a distance modulus 
of Vo ~ Mv = 9.3±0.5 (725 pc) and a metallicity of [Fe/H] = -0.24±0.27. Through isochrone fitting we find 
an age of logi = 8.6±0.1 (400 Myr). Our optical photometry and JHK from 2MASS are combined to derive 
effective temperatures of cluster member stars. The effective temperature distribution along the main sequence 
of the cluster shows several gaps. A test to study the significance of these gaps in the main sequence of the HR 
diagram has been applied. The method is also applied to several other open clusters (Pleiades, Hyades, NGC 1817 
and M 67) to construct a sequence of metallicities and ages. The comparison of the results of each cluster gives 
four gaps with high significance (one of them, centred at 4900 K, has not been previously reported). 
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1. Introduction photometry in the Sloan system with a magnitude limit 

g'~l8. A recent study by Wu et al. <|20Tj5|> gives BATC 

The open cluster NGC 2548 (C0811-056), also known as photometry in 13 filters . Radi al velocities have been stud- 

M 48, in Hydra (a 200 o=8 h 13 m 48 s , 5 200 o =__ 5 ° 48 ') with ied by Geyer & Nelles lfl» with a focal reduced spec- 

an esti mated distance of 630 pc (Pesch [I96IJ or 530 pc trograph. They give data for 23 stars but with a very low 

(Claria[I985J) , is aninter mediate- age open cluster, around qualit y The only quality measurement is from Wallerstein 

logi = 8.5 (Lynga[1987J), with a slightly poorer CN abun- et al fit%ft f star WEBDA 1560. 
dance than the giants of the Hyades but significantly richer 
than the K giants of the solar neighborhood (Claria 1985 
Twarog et al. I1997f) . It has been very poorly studied in 



This study of NGC 2548 is part of a series of astromet- 
ric and photometric analyses of open clusters that follows 

. , . , , , • • , , already published results on NGC 1817 (Balaguer-Nuhez 

spite of being an extended obiect with an apparent diame- . . ; W ,,,, ■ ■ . , . .... . . , .,, 

J~ ~f Qrv rrC * wi-ram* ™ * *a> ^n^^iram et aL MM [2004b| and that will be completed with a 



study of M 67 (Balaguer-Nuhez et al. 



ter of 30' (Trumpler lTT%U|l or even 54' (Collinder lllEfl)! and 
brilliant enough to be in the Messier list (XVIII century) 
as number 48 (Messier lTMTI) . It was even considered inex- Absolute proper motions of 501 stars within a 1?6 x 
istent for several years due to the fact that Messier quoted 1?6 area in the NGC 2548 re g ion > from automatic MAMA 
its coordinates with an error of 5°. Rider et al. itSiTOI) gives measurements of 10 plates taken with the double astro- 
graph at Z6-Se station of Shanghai Observatory, were 
studied by our group (Wu et al. 2002, hereafter Paper I). 



Send offprint requests to: Balaguer-Nufiez, L., e-mai 

Lola.Balaguer@am.ub.es In this paper we discuss the results of a deep CCD pho- 

* Tables 3, 4 and 5 are only available in electronic form at the tometric study of NGC 2548, covering an area of 34' x 34' 

CDS viaanonymousftptocdsarc.u-strasbg.fr (130.79.128.5) or down to V f» 22. Section 2 contains the details of the 

via http://cdsweb.u-strasbg.fr/cgi-bin/qcat?J/A+A/ CCD observations and their reduction and transforma- 
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tion to the standard system. In Sect. 3 we discuss a new 
astrometric membership segregation based on the compar- 
ison between parametric and non-parametric approaches 
applied to the proper motions from Paper I. In Sect. 4 
we discuss the colour-magnitude diagram and identify the 
sample of candidate cluster members using astrometric as 
well as photometric criteria. Section 5 contains the deriva- 
tion of the fundamental cluster parameters of reddening, 
distance, metallicity and age. In Sect. 6 we calculate ef- 
fective temperatures and study the significance of gaps 
in the main sequence. The method is also applied to our 
results on NGC 1817 and NGC 2682 (M 67) and to the 
well-studied Hyades and Pleiades clusters, allowing us to 
test the presence of the gaps as a function of age and 
metallicity. Section 7 summarizes our conclusions. 

2. The data 

2.1. Observations 

Deep Stromgren CCD photometry of the area was per- 
formed at the Calar Alto Observatory (Almeria, Spain) in 
January 1999 and January 2000 using the 1.23 m telescope 
of the Centro Astronomico Hispano-Aleman (CAHA) and 
in January 1999 and February 2000 using the 1.52 m tele- 
scope of the Observatorio Astronomico Nacional (OAN). 
Further data were obtained at the Observatorio del Roque 
de los Muchachos (ORM, La Palma, Canary Islands, 
Spain) in February 2000 using the 2.5 m Isaac Newton 
Telescope (INT) of ING (equipped with the Wide-Field 
Camera, WFC), and in December 1998 and February 2000 
using the 1 m Jakobus Kapteyn Telescope (JKT) of ING, 
with the Hp filter. The poor quality of the images ob- 
tained in the 1998/99 runs and in the OAN 2000 observa- 
tions, due to adverse meteorological conditions, prevented 
us from making use of the data collected during those 
nights. A log of the observations with the total number of 
frames, exposure times and seeing conditions is given in 
Table H] 

We obtained photometry for a total of 4806 stars in 
an area of 34'x34' around NGC 2548, down to a limiting 
magnitude V~22. The area covered by the observations is 
shown in the finding chart of the cluster (Fig.^). Due to 
the lack of Hp filter at the WFC-INT, it was only possible 
to measure it at the JKT and CAHA telescopes, thus lim- 
iting the spatial coverage with this filter. Only 283 stars 
have Hp values. Of those only 253 have values of the rest 
of the Stromgren filters. 

Beside long, deep exposures, additional shorter expo- 
sures were obtained in order to avoid saturation of the 
brightest stars. 

2.2. Reduction and transformation to standard system 

The reduction of the photometry is explained at length 
in Balaguer-Nuhez et al. ( 2004b). Our general procedure 
has been to routinely obtain twilight sky flats for all the 
filters and a sizeable sample of bias frames (around 10) be- 



Fig. 1. Finding chart of the area under study. The covered 
area is marked in black on an image of a plate plotted with 
Aladin. 

fore and/or after every run. Flat fields are typically fewer 
in number, from five to ten per filter. Two or three dark 
frames of 2000 s were also taken. IRAF 1 routines were 
used for the reduction process. The bias level was evalu- 
ated individually for each frame by averaging the counts 
of the most stable pixels in the overscan areas. The 2-D 
structure of the bias current was evaluated from the aver- 
age of a number of dark frames with zero exposure time. 
Dark current was found to be negligible in all the cases. 
Flatfielding was performed using sigma clipped, median 
stacked, dithered twilight flats. 

Our fields are not crowded. Thus, synthetic aperture 
techniques provide the most efficient measurements of rel- 
ative fluxes within the frames and from frame to frame. 
We used the appropiate IRAF packages, and DAOPHOT 
and D AO GROW algorithms fStetson lil%7llTMi|l . We an- 
alyzed the magnitude growth curves and determined the 
aperture correction with the IRAF routine MKAPFILE. 

For the WFC images from the INT, we employed 
the pipeline specifically developed by the Cambridge 
Astronomical Survey Unit. The process bias subtracts, 
gain corrects and flatfields the images. Catalogues are 
generated using algorithms described in Irwin (1985). 
The pipeline gives accurate positions in right ascen- 
sion and declination linked to the USNO-2 Catalogue 
(Monet et al. 1998), and instrumental magnitudes 
with their corresponding errors. A complete descrip- 

1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation. 
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Telescope 


Date 


Seeing(") 


n. of frames 




Exp. 


Times 


(sec) 












u 


V 


b 


y 


Hp 


1.23 m CAHA 


1999/01/12-15 


(1) 


13 


1900 


800 


400 


400 


2000 


1.23 m CAHA 


2000/01/05-10 


1.1 


21 


2200 


1400 


900 


800 


1400 


1.52 m OAN 


1999/01/13-16 


(1) 


15 


1900 


800 


400 


400 


2000 


1.52 m OAN 


2000/02/07-14 


(1) 


5 










1400 


1 m JKT 


1998/12/11-14 


(1) 


26 


2000 


1200 


800 


700 


1200 


1 m JKT 


2000/02/02-06 


1.1 


18 










2000 


2.5 m WFC-INT 


2000/02/02-03 


1.3 


17 


2000 


2000 


1200 


500 





(1) Poor weather conditions. Images not used in the final data. 

tion can be found in Irwin & Lewis (2001) and in 
http : //www . ast . cam. ac . uk/~wf c sur/ index .php 

The coefficients of the instrumental-to standard trans- 
formation equations were computed by a least squares 
method using the instrumental magnitudes of the stan- 
dard stars and their standard magnitudes and colours in 
the uvby — Hp system. Up to 68 standard stars from the 
cluster M 67 (Nissen et al. n^57|l were observed depend- 
ing on the size of the field. Four to six short exposures 
in every filter were taken every night with a magnitude 
limit of V—18. Those standard stars with residuals greater 
than 2a were rejected. Typically that involves a 10-15% of 
the standard stars, mainly variables from the M 67 field. 
The reduction was performed for each night independently 
and in two steps. The first step is to determine the extinc- 
tion coefficients for each passband from the standard stars. 
With the extinction coefficients fixed, the transformation 
coefficients to the standard system were fitted. 

The mean errors as a function of apparent visual mag- 
nitude are given in Table for the NGC 2548 stars. 

Table 3 lists the u, v, b, y, Hp data for all 4806 stars 
in a region of 34' x 34' around the open cluster NGC 2548 
(Fig. ^). Star centres are given in the frame (x,y) and 
equatorial (aj2ooo A2000) coordinates. An identification 
number was assigned to each star following the order of 
increasing right ascension. Column 1 is the ordinal star 
number; columns 2 and 3 are aj2ooo and <5j2ooo; columns 
4 and 5 are the respective x, y coordinates in arcmin; 
columns 6 and 7 are the (b — y) and its error, 8 and 9 the 
V magnitude and its error, 10 and 11 the mi and its error, 
12 and 13 the c\ and its error, and 14 and 15 the Hp and 
its error. In column 16, stars considered candidate mem- 
bers (see Sect. 4.1.) are labelled 'M', while those classified 
as non- members show the label 'NM'. 

The cross-identification of stars in com- 
mon with the astrometry (Paper I), WEBDA 
(http://obswww.unige.ch/WEBDA), Hipparcos (ESA, 
ITM7|>, T ycho-2 (H0g et al. l2tjOU|l and USNO-2 (Monet et 
al. I1998fl catalogues is provided in Table 4. 



2.3. Comparison with previous photometry 

Only three stars in the NGC 2548 area have been pre- 
viously studied using Stromgren photometry (WEBDA 
0366, 1560, 2156). Unfortunately, none of them is inside 
the area covered by our photometry. Only one of them 
(WEBDA 1560) is considered a cluster member from our 
astrometry and will be studied in Sect. 5. 

The V magnitude derived from the y filter can be com- 
pared with the published broadband data. There are 30 
stars in common with Pcsch (1961). The corresponding 
mean difference in V, in the sense ours minus Pesch' is 
-0.01±0.03. 

Transformations between B — V and b — y from sev- 
eral authors (see Moro & Munari, I2()()()|l fails to cover 
the whole range under study. We can find a linear rela- 
tion between the two indices: B — V = (1.632±0.031)(6 — 
y) - (0.038 ± 0.010), N=30. The standard deviation of 
the residuals about the mean relation is 0.048, where the 
typical uncertainty in B — V is 0.02 and in b — y is 0.016. 

3. Astrometric analysis 

Paper I gives absolute proper motions of 501 stars within 
a 1?6 x 1?6 area in the NGC 2548 region up to B inst ^U, 
from automatic MAMA measurements of 10 plates and a 
maximum epoch difference of 82 years. Membership de- 
termination is calculated in Paper I with a 9-parametric 
Gaussian model and a list of stars with a probability 
higher than 0.7 gives 165 cluster members. 

NGC 2548 is a very extended object with a complex 
structure with a double core, prolate shape and a tidal 
tail with a clump (Bergond et al. I2001J1 . It has been sug- 
gested that this secondary clump is associated with the 
last strong disk shock that occurred between 20 and 40 
Myr ago. Confirmation of members at large radii will trace 
the distribution of stars that are currently leaving the clus- 
ter. This will help to constrain models of the tidal disrup- 
tion of open clusters. 

To complement the cluster/field segregation analysis 
of the astrometry from Paper I we have applied a non- 
parametric method to the proper motion data, as ex- 
plained in Balaguer-Nunez et al. f 2"004aj ■ In the non- 
parametric approach we are able to differentiate the clus- 
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Table 2. Number of stars observed (JV) and mean internal errors {a) as a function of apparent visual magnitude. 
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8- 9 


3 


0.019 


3 


0.028 


2 


0.016 


2 


0.002 


2 


0.017 


9-10 


19 


0.012 


19 


0.011 


19 


0.013 


16 


0.019 


5 


0.003 


10-11 


25 


0.014 


25 


0.016 


25 


0.016 


24 


0.010 


6 


0.009 


11-12 


32 


0.012 


31 


0.013 


31 


0.018 


31 


0.016 


10 


0.030 


12-13 


45 


0.012 


43 


0.013 


43 


0.016 


43 


0.011 


12 


0.017 


13-14 


71 


0.012 


71 


0.015 


69 


0.020 


69 


0.017 


17 


0.027 


14-15 


152 


0.012 


152 


0.023 


146 


0.029 


146 


0.030 


29 


0.016 


15-16 


260 


0.015 


260 


0.034 


259 


0.047 


251 


0.047 


42 


0.025 


16-17 


380 


0.022 


380 


0.043 


376 


0.061 


360 


0.063 


53 


0.033 


17-18 


514 


0.030 


514 


0.035 


501 


0.043 


468 


0.043 


51 


0.036 


18-19 


639 


0.027 


639 


0.028 


592 


0.031 


477 


0.032 


8 


0.046 


19-20 


752 


0.012 


752 


0.015 


702 


0.022 


383 


0.038 






20-21 


859 


0.023 


859 


0.028 


662 


0.048 


160 


0.052 






21-22 


763 


0.049 


763 


0.063 


262 


0.089 


21 


0.082 






22-23 


222 


0.095 


222 


0.122 


28 


0.139 










Total 


4736 




4733 




3717 




2451 




235 





ter population without the need for any a priori knowl- 
edge. Following Galadf-Enrfquez et al. (1998), we can per- 
form an empirical determination of the probability den- 
sity functions (PDFs) evaluating the observed local den- 
sity in each node of a two-dimensional grid in the vector 
point diagram (VPD), with a normal circular kernel. The 
smoothing parameter h (Gaussian dispersion) is chosen 
using Silverman's rule H198fi[) . The procedure was tested 
for several subsamples applying different proper motion 
cutoffs. Satisfactory results are obtained with a proper 
motion cutoff of < 15 mas yr _1 . 

The empirical frequency function determined from the 
VPD corresponding to the area occupied by the cluster 
is made up of two contributions: cluster and field. To dif- 
ferentiate the two populations we need to estimate the 
field distribution. For this purpose, we studied the VPD 
for the plate area outside a circle centered on the cluster. 
The centre of the cluster was chosen as the point of highest 
spatial density. We did tests with circles of very different 
radii, searching a reasonable tradeoff between cleanness 
(absence of a significant number of cluster members) and 
signal-to-noise ratio (working area not too small) . The ker- 
nel density estimator was applied in the VPD to these 
data, yielding the empirical frequency function, for a grid 
with cell size of 0.2 mas yr _1 , well below the proper mo- 
tion errors. 

We finally found that the area outside a circle with 
a radius of 35' centered on the cluster yields a clean fre- 
quency function with low cluster contamination and low 
noise. The empirical frequency functions can be normal- 
ized to yield the empirical PDFs for the mixed population 
(circle) , for the field (outside the circle) and for the cluster 
(non-field) population. Figure displays these three func- 
tions. Of course, the field PDF estimated in the outer area 
cannot be an absolutely perfect representation of the true 
field PDF in the whole area. This introduces undesired 



noise in the frequency function of the cluster, as shown 
by the negative density values found in several zones and 
by some positive fluctuations without physical meaning. 
These negatives values allow us to estimate the typical 
noise level, 7, present in the result. To avoid meaning- 
less probabilities in zones of low density, we restricted 
the probability calculations to the stars with cluster PDF 
> 37. The maximum of the cluster PDF is located at 
(fj, a cos 5, (X S ) = (-1.2±0.2,2.2±0.2) mas yr" 1 . 

However, this cluster proper motion is different to the 
value found in Paper I. As discussed by many authors 
(Galadf-Enrfquez et al. 119981 for example) one of the lim- 
itations of the parametric approach is the trend of the 
circular Gaussian distribution, used to fit the cluster, to 
assume an excessive width to improve the representation 
of the field distribution. The cluster mean proper motion 
will then be thus affected. To measure the influence of this 
effect, we decided to apply the parametric method but to 
fix the internal velocity dispersion of the cluster at zero. 
This way the model will assign to the cluster Gaussian dis- 
tribution a width related only to measurement errors. We 
then obtained a mean proper motion of (/i a cos<5, us) = 
(-1.10±0.08,2.09±0.08) mas yr" 1 , in agreement with the 
cluster proper motion obtained by the non-parametric ap- 
proach. 

The non-parametric technique does not take into ac- 
count the errors of the individual proper motions, there- 
fore it does not make any particular distinction be- 
tween bright or faint stars, different epoch spread and 
so on. However, the FWHM of the empirical cluster 
PDF provides an estimation of the errors of the distri- 
bution. We obtained a FWHM of ~4.2±0.2 mas yr" 1 . 
If the Gaussian dispersion owing to the smoothing pa- 
rameter h — 1.44 mas yr _1 is taken into account, this 
FWHM corresponds to a value of 1.53 mas yr _1 . But from 
Paper I we know that the mean proper motion precision 
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is 1.18 mas yr _1 which gives us an intrinsic dispersion 
component of 0.97 mas yr _1 , (3 km s _1 at the distance of 
725 pc from Sect. 5), of the same order but slightly lower 
than the value obtained by the membership determina- 
tion in Paper I. This indicates that the intrinsic velocity 
dispersion of the cluster cannot be neglected. Although 
fixing it to zero improves the determination of the mean 
motion of the cluster (position of the centre of the fitted 
Gaussian), the membership probability results are more 
meaningful taking into account the intrinsic dispersion. 
In our analysis, we will use the parametric results from 
Paper I. However, slight differences in the center of the 
adopted Gaussian do not affect the segregation, since the 
stars with highest probability of being members are almost 
the same in both cases. 

The cluster membership probability histogram (Fig.|3J) 
shows a clear separation between cluster members and 
field stars in both approaches: the solid line is the tradi- 
tional parametric method (from Paper I) while the dotted 
line is the non-parametric approach. The non-parametric 
approach to cluster/field astrometric segregation gives us 
an expected number of cluster members from the inte- 
grated volume of the cluster frequency function in the 
VPD areas of high cluster density, where PDF> 37. This 
integration predicts that the sample contains 91 cluster 
members. Sorting the sample in order of decreasing non- 
parametric membership probability, P/yp, the first 91 stars 
are the most probable cluster members. The minimum 
value of the non-parametric probability (for the 91-st star) 
is P NP = 82%. Table 5 lists the P NP for the 501 stars. 

There is not an equivalent rigorous way to decide where 
to set the limit among members and non-members in the 
list sorted in order of decreasing parametric membership 
probability, Pp. But, if we accept the size of the cluster 
predicted by the non-parametric method, 91 stars, we can 
consider that the 91 stars of highest Pp are the most prob- 
able members, according to the results of the parametric 
technique. The minimum value of the parametric proba- 
bility from Paper I (for the 91-st star) is Pp — 92%. 

With these limiting probabilities (Pnp > 0.82; Pp > 
0.92), we get a 91% (458 stars) agreement in the segre- 
gation yield by the two methods. The 43 remaining stars 
(9%) with contradicting segregation should be carefully 
studied. Discrepancies among the two approaches are ac- 
tually expected due to the statistical nature of the meth- 
ods themselves. 

As in Balaguer-Nunez et al. (2004a), to set up a final 
and unique list, and trying not to reject true members, 
we accept as probable members of this cluster those stars 
classified as members by at least one of the two meth- 
ods. This is equivalent to merging both lists (each with 91 
stars) and eliminating duplicated entries. This way we get 
a list of 118 probable member stars. 

As commented on the Introduction, the available in- 
formation on radial velocities is not accurate nor complete 
enough to be useful in improving the membership segre- 
gation. 



We have studied twenty stars in the area of the sec- 
ondary clump on the tidal tail proposed by Bergond et al. 
(2001). Only three of them appear to be cluster members. 
The rest are randomly distributed in the VPD. Because 
our limiting magnitude for astrometry _B; ns t«14 is brighter 
than the B cu t«14.8 of Bergond et al. (2QCJTJ), we cannot 
discard the reality of this clump. Unfortunately, our pho- 
tometry does not cover this area. Deeper studies will be 
necessary. 

4. Colour-Magnitude diagrams 

We use the Vvs (v — y) colour-magnitude diagram for our 
study. The colour-magnitude diagram based on this colour 
index defines the main sequence of a cluster significantly 
better than the traditional Vvs (b — y) diagram (Fig.0]lcft 
and centre) . The colour-magnitude diagram of all the stars 
in the area displays a fairly well-defined main sequence. 

4.1. Selection of candidate member stars 

Our photometric measurements help to reduce the possi- 
ble field contamination in the proper motion membership 
among bright stars, as well as to enlarge the selection of 
members towards faint magnitudes. Among those astro- 
metric member stars with photometric measurements, we 
find one star that is not compatible with the sequence 
of the cluster outlined in the colour-magnitude diagram. 
Our astrometric segregation of member stars has a limit- 
ing magnitude of V ~ 13. From this point down to V = 18 
we construct a ridge line following a fitting of the obser- 
vational ZAMS (Cra wford UTTTBI H"5?8l H3T91 Hilditch et 
al. 113551 Olsen ll^Mjl in the V — (v — y) diagram. A se- 
lection of stars based on the distance to this ridge line is 
then obtained. The chosen margin for candidates includes 
all the stars between V + 0.5 and V — 1 from the ridge 
line, as shown in the right panel of Fig. The margins 
were chosen to account for observational errors and the 
presence of multiple stars. 

This preliminary photometric selection is refined in the 
colour-colour diagrams (Fig.EJ with the help of the stan- 
dard relations from the same authors. A final selection of 
331 stars in the area is plotted in Fig. as empty circles 
in the [mi] - [c x ], mx - (b - y), c x - (b - y) and [c\] - Hp 
diagrams. 

5. Physical parameters of the cluster 

The stars selected as possible cluster members were clas- 
sified into photometric regions and their physical param- 
eters were determined as explained in Balaguer-Nuhez et 
al. <|2004b|) . The algorithm, described in Masana et al. 
(2005a) and Jordi et al. (1997), uses uvby — Hp photom- 
etry and standard relations among colour indices for each 
of the photometric regions of the HR diagram. Absolute 
magnitude, effective temperature and gravity as well as 
the corresponding reddening, distance modulus, metallic- 
ity and raw spectral type and luminosity class are calcu- 
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Fig. 4. The colour-magnitude diagram of the NGC 2548 area. Empty circles in the right figure are the astrometric 
members. Thick line is a shifted ZAMS, with the chosen margin for candidate members (V + Q.5,V — 1) in thin lines. 
See text for details. 



lated for each star. Typical errors are 0.25 mag in My, 
0.15 dex in [Fe/H], 270 K in T cff , 0.18 dex in logg and 
0.015 mag in E(b — y). 

Only 62 stars among the 331 candidate members have 
Hp measurements. Thus, the computation of physical pa- 
rameters is only possible for that subset. The results are 
shown in Fig. Excluding peculiar stars and those with 
inconsistencies in their photometric indices and applying 
an average with a 2a clipping to that subset, 39 stars re- 
main. We found a reddening value of E(b — y) — 0.06±0.03 
(corresponding to E(B — V) = 0.08) and a distance mod- 
ulus of V - My = 9.3±0.5 (725 pc, i.e. about 200 pc 
above the galactic plane). Metallicity is better calcu- 
lated studying only the 26 F and G stars in our sam- 
ple following Masana et al. I|2005aj) . We found a value of 
[Fe/H] = -0.24±0.27. 

We have examined the bibliography in search of bright 
stars in the area of NGC 2548 and have found only 
one compatible with being a member of the cluster: The 
brightest star (WEB DA 1560), that was saturated in our 
data. We therefore took its values from Olsen ( 1993). The 
stars in the red giant region are detailed in Tabled Three 
of them are known spectroscopic binaries. We are cautious 
about our values of star WEBDA 1218, as the known val- 
ues of V are 9.64 fPesch lTMTf and 9.63 (Oia lT9T6) . about 
0.2 mag brighter. Star WEBDA 0870 was considered to 
be an astrometric non-member by Ebbinghausen (1939) 
and in Paper I, while Claria (1985) studied its weak CN 
peculiarity and opened the possibility of a cluster mem- 
ber with peculiar composition and proper motion. Star 
WEBDA 1241 was considered to be a cluster member in 
Paper I. 

The recent publication by Clem et al. ( 2003 °f ern ~ 
pirically constrained colour-temperature relations in the 
Stromgren system enables an isochrone fitting to our re- 
sults. The best fit is found for the Pietrinferni et al. 
I|2004|) tracks. Figure0shows isochrones of [Fe/H] = —0.25 



for canonical models (right panel) and for models with 
overshooting (left panel). Taking into account the num- 
ber of stars before and after the MS hook, we found a 
best fit for models with overshooting. The estimated age 
is of 400±100 Myr (logt = 8.6). Using a different set 
of tracks by Schaller et al. (|1992f) we also found agree- 
ment with a best estimation of the age of logt = 8.6±0.1. 
Neither Pietrinferni et al. I|2004j> nor Schaller et al. (|1992|l 
isochrones provide a perfect fit to the giant members. 
Similar discrepancies were found by Lapasset et al. (2000) 
studying the intermediate-age cluster NGC 2539. 

Our results are consistent with previous studies. Pesch 
([13611) gives an E(B - V) = 0.04±0.05 and a distance of 
630 pc from U BV photoelectric observations of 30 stars. 
Claria (1985) studied DDO photometry of 5 giant stars 
and gives a value of E(B — V) = 0.06 and a distance of 
530 pc. Twarog et al. Ijl997|l revised those values and gave 
an E{B-V) =0.05, a distance of 590 pc and [Fe/H] =0.08 
from 3 stars. Harris l|1976fl gives a value of logt = 8.28 
from an analysis of nine stars. H0g & Flynn (1998) give 
[Fe/H] = -0.13, E(B - V) = 0.04 and logt = 8.50 from 
DDO photometry of three giants. Rider et al. :|200 1 ■ fix a 
value E(B — V) = 0.03 from bibliography and give values 
of distance 700 pc, [Z/Zq] = 0.0 and an age of 400 Myr 
by isochrone fitting to u'g'r'i'z' photometry. The recent 
study by Wu et al. (2005 ) gives a slightly larger distance 
(780 pc) and a lower age (0.32 Gyr) also fitting isochrones 
with E(B — V) = 0.04 and solar metallicity. A summary 
of this information is given in Table for clarity. 

6. Gaps in the Main Sequence 

Gaps in the main sequence have been observed and de- 
bated for a long time. In the Hyades, de Bruijne et al. 
(2000) confirmed the reality of two gaps near B — V ~ 0.38 
and B - V ~ 0.48 (T cff ~ 6400 K). Evidence for the ex- 
istence of gaps among field stars has also been debated 



Balaguer- Nunez et al: uvby — Hp Photometry and Membership of NGC 2548 



7 



Table 6. Red giants in NGC 2548, with their membership. The first five are from the WEBDA list of Red Giants. 
Three of them are known spectroscopic binaries. 



Phot. Id. 


WEBDA 


(& — y) 


V 


Pp 


Pnp 


M/NM 




1006 


0870 


0.670±0.0185 


9.712±0.0185 


0.07 


0.84 


M? 




2556 


1218 


0.335±0.0370 


9.862±0.0305 


0.83 


0.85 


M 




2732 


1260 


0.370±0.0002 


9.050±0.0047 


0.97 


0.84 


M 


SB 


2894 


1296 


0.504±0.0004 


9.215±0.0023 


0.94 


0.86 


M 


SB 


(1) 


1560 


0.724±0.0020 


8.199±0.0060 


0.88 


0.84 


M 


SB 


3971 


1521 


0.312±0.0079 


9.671±0.0081 


0.00 


0.00 


NM 




4557 


1616 


0.719±0.0045 


9.890±0.0045 


0.00 


0.00 


NM 




2645 


1241 


0.777±0.0072 


9.565±0.0072 


0.81 


0.66 


NM? 





(1) Photometry taken from Olsen (1993) 




Fig. 7. Isochrones from Pietrinferni et al. (2004) for scaled solar models of [Fe/H] = —0.25 and ages of 300, 400, 500, 
600 Myr for models with overshooting (left panel) and 200, 300, 400 Myr for canonical models (right panel). Empty 
circles are candidate members. Red giants are analyzed in detail (see text): diamonds are spectroscopic binaries, squares 
are astrometric non-members. The adopted reddening and distance modulus are E(b — y) = 0.06 and Vq — My = 9.3. 



(Newberg & Yannv IT998f) . and recently Kovtyukh et al. 
( 2004) showed a new gap among field stars in the effective 
temperature range 5560 to 5610 K. 

6.1. Method 

The characteristics of our photometry and member selec- 
tion procedure allow us to pursue the search for gaps in 
the main sequence not only for the previously discussed 
cluster, NGC 2548 (M 48), but also in the other clusters 
included in our project: NGC 1817 (Balaguer-Nunez et al. 
I20fl4al I2004b|l and NGC 2682 (M 67, Balaguer-Nunez et 
al. E005). To check the reliability of the results, we ap- 
ply the same gap search to the Pleiades (M 45) and the 
Hyades, two well-known clusters in which gaps have been 



previously reported (de Bruijne et al. 12)00, Belikov et al. 
119981 and specially their Figs. 1 and 3, respectively.) 

We estimate effective temperatures for our selection of 
candidate members following a new approach (Ribas et al. 
I2TKH Masana et al. I2005bjl based on fitting observed IR 
photometry with accurately calibrated synthetic photom- 
etry. We use our values of V magnitude and the 2MASS 
values for J, H , and K magnitudes. The process requires 
metallicity and surface gravity data that we compute from 
our Stromgren photometry as explained in previous sec- 
tions. The method is restricted to the temperature interval 
from 4000 K to 8000 K, with quoted uncertainties of 0.5- 
1.3%. The upper temperature limit is due to the increased 
dependence of the results on the accuracy of logy, while 
the lower limit is related to the decreasing performance of 
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Table 7. Comparison of the physical parameters of the cluster with results from other authors. See the text for 
detailed information of the assumptions made in each case. 



Reference 


E(B - V) 


[Fe/H] 


log* 


Distance 


Comments 




mag 


dex 




pc 




This work 


0.08 


-0.24±0.27 


8.6±0.1 


725 




Wu et al. (1217051 


0.04 


0.0 


8.51 


780 


Isochrone fit 


Rider et al. 150011 




0.0 


8.60 


700 


Isochrone fit 


H0g & Flynn (THJQgt 


0.04 


-0.13 


8.50 




DDO Phot. 3 stars 


Harris 






8.28 




Phot. 9 stars 


Twarog et al. itTOOTl 


0.05 


0.08 




590 


DDO Phot. 3 stars 


Claria 115551 


0.06 


0.1 




530 


DDO Phot. 5 stars 


Pesch itMflTl 


0.04 






630 


UBV Photoelectric 30 stars 



the models because of molecular bands. Ribas et al. ( 2003 ) 
show that, in this range, the procedure is essentially insen- 
sitive to the adopted value of [Fe/H] and logg: uncertain- 
ties of 0.3 dex in metallicity and 0.5 dex in gravity induce 
T e ff deviations inferior to 0.5%. The number of member 
stars from our selection of clusters that also have 2MASS 
photometry in the range of temperatures under study is of 
269 stars for NGC 2548, 307 stars for NGC 1817 and 588 
stars for NGC 2682. For Pleiades we obtain 225 member 
stars and for Hyades 76. 

After getting the T e g values, we proceeded to the gap 
search following a method analogous to that proposed by 
Rachford & Canterna JgDDUf). A simple \ 2 test with one 
degree of freedom is used to evaluate the significance of 
any candidate gap. To do this, we take a candidate gap of 
width Wi n and we compute the number of stars within it, 
-/V; n . This number is compared with the stars located on 
both sides of the gap. We take two bins of equal width, 
Woutj at the sides of the candidate gap, and we count the 
stars inside them, rt ou t. Then, n out is scaled to the size of 
the gap to give N out = Wi n n out /(2Wo U t) and this quantity 
is compared to iVj n to get the x 2 value. To limit the effect 
of small number statistics, and again following Rachford 
& Canterna (2000), the computation is done only when at 
least five stars are present on each side of the gap, which 
guarantees that n out > 10. 

We performed this calculation for a grid of gap widths 
75 K < W m < 500 K, at intervals of 25 K, and placing the 
gap centres every 1 K over the intersection of the tempera- 
ture range covered by the method with the range covered 
by the photometry. To prevent edge effects, intervals of 
300 K were avoided at the extremes of the intersection. 

After several trials, we chose W ou t — 100 K for our 
photometry, but the Pleiades and Hyades required W out — 
200 K, due to the smaller size of the photometric sample. 
It is important to note that widening W ou t can alter the 
relative value of x 2 and its associated probability, but with 
little effect on the significance of the gap centre position 
or width. 

Having computed x 2 f° r the whole range of gap centres 
and sizes, we selected the local maxima of the resulting 
table. We observed that changes in W n do not make gaps 



appear or disappear, but only affect their significance. All 
significant gaps have probability values higher than 0.99. 

6.2. Results 

Table IS] gives the centres and widths (Win) for the signifi- 
cant gaps in the temperature range studied. The two well- 
known Bohm-Vitense gaps are clearly detected in Hyades 
and also Pleiades, which serves as a check of the reliabil- 
ity of the method. Both gaps are significant in the whole 
sample of clusters analysed, with the exception of M 67, 
due to the much higher age of this stellar system (the 
main sequence is so evolved that the corresponding effec- 
tive temperature range is no longer populated). 

The gap recently reported by Kovtyukh et al. (2004) 
in field stars at 5560-5610 K, and already suspected in the 
Hyades (de Bruijne et al. .2000), stands out not only in 
the Hyades, but also in the other clusters surveyed. In the 
case of M 67, Fan et al. (1992J) also raised the possible 
existence of this gap placed around 1 M . 

The reliability of the new gap detected around T c ff ~ 
4900 K is stressed by its appereance in all the clusters 
under study with the only exception of NGC 1817. The 
reason is that, in this case, our photometry does not reach 
stars colder than 5000 K. 

Even though we know that some amount of field con- 
tamination is present in our candidate member selection, 
it is worth noting that the significance of the four gaps in 
our photometric data is outstanding. If an even more reli- 
able member selection would be possible, this could make 
the gaps even more evident in our data. 

The position and width of the gaps are compared 
with the age and metallitity of the clusters in Fig. [SJ No 
clear trend can be drawn from this comparison. However, 
metallicity effects could be masked due to the uncertain- 
ties in the metallicity determinations used for NGC 2548, 
NGC 1817 and NGC 2682 (M 67) that are of the order of 
0.2 dex. 

Independently of the physical explanation for the ex- 
istence of these gaps, there is no reason, in principle, to 
expect them to appear at the same positions and with 
similar sizes in all the clusters. There can be a com- 
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Cluster 


Age 


[Fe/H] 


Centre 


Width 


Centre 


Width 


Centre 


Width 


Centre 


Width 




Gyr 


dex 


K 


K 


K 


K 


K 


K 


K 


K 


Pleiades (M 45) 


0.1 


-0.03 


7377 


300 


6821 


325 


5577 


275 


4845 


275 


NGC 2548 (M 48) 


0.4 


-0.24 


7099 


475 


6305 


425 


5465 


350 


4785 


350 


Hyades 


0.7 


0.15 


7006 


475 


6427 


250 


5452 


225 


4972 


275 


NGC 1817 


1.1 


-0.34 


7323 


175 


6674 


350 


5701 


300 






NGC 2682 (M 67) 


4.2 


0.01 










5593 


125 


5017 


200 



Table 9. Oricntative and approximate characterization of 
the mean loci of the gaps detected. 



< T off > (K) 


ff<T rff > (K) 


B-V 


Spec. 


Mass (M ) 


7200 


180 


~0.30 


F0 


~1.6 


6600 


230 


~0.40 


F2-F5 


~1.3 


5600 


100 


~0.70 


G5-G8 


~0.80 


4900 


110 


~0.90 


K2 


~0.75 



plex dependency on metallicity, age and other parame- 
ters. However, the clusters studied show some regularity 
that allows to distinguish, as mentioned, four independent 
gaps. Even though we know that this kind of average may 
lack physical meaning, we give a table with mean, approx- 
imate parameters that caracterize the loci of the four gaps 
(Table 0). The colour indexes and masses are taken from 
standard relations (Schmidt-Kaler 1982| and assume solar 
metallicity. 

Several theoretical explanations have been proposed 
for the three hotter gaps (see the references already given). 
Most but not all, point to the need to improve the 
treatment of rotation and convection in stellar models. 
However, theoretical and/or empirical references to the 
fourth, colder gap are lacking. Spectroscopy could confirm 
or reject the reality of this gap whose significance is simi- 
lar to that of the others, on the basis of our photometric 
data. 

7. Conclusions 

In this paper we give a catalogue of accurate Stromgren 
uvby — Hp and J2000 coordinates for 4806 stars in an area 
of 34'x34' around NGC 2548. 

We give a selection of probable members of NGC 2548, 
combining the photometric study with an astrometric 
analysis using parametric as well as non-parametric ap- 
proaches. An improved determination of the physical pa- 
rameters of this cluster based on our photometry gives: 
E(b -y) = 0.06±0.03, [Fe/H] = -0.24±0.27, a distance 
modulus of Vq — My — 9.3±0.5 and an age of logi = 
8.6±0.1. The values are consistent with previous studies 
(see Chen et al. 2003 and references therein). 

We perform a search for possible gaps in the main se- 
quence of this cluster and, also, of the other clusters anal- 
ysed in our project (NGC 1817, M 67). The Pleiades and 
Hyades are added as a check and to construct a sequence 



of ages and metallicities. We find four gaps in the effec- 
tive temperature interval 4000-8000 K. The coldest gap, 
around 4900 K, had not been previously reported. 
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Fig. 2. Empirical probability density functions in the 
kinematic plane. Top: mixed sample from the inner circle 
of 35'. Center: field population from outside this circle. 
Bottom: cluster population of NGC 2548 
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Fig. 3. The histogram of cluster membership probability 
of NGC 2548. The solid line gives the results for tradi- 
tional parametric method (Paper I), while the dotted line 
corresponds to the non-parametric approach. The arrows 
mark the limiting probabilities for member selection for 
each method. 




Fig. 5. The colour-colour diagrams of NGC 2548. Empty 
circles denote candidate members of NGC 2548, chosen 
with astrometric and non-astrometric criteria as explained 
in Sect. 4.1. The thick line is the standard relation shifted 
by E(b — y) = 0.06 when necessary. 
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Fig. 6. The histograms of the distance modulus, red- 
dening and metallicity of the selected member stars of 
NGC 2548 with Hp measurements. The arrows indicate 
the mean values adopted for the cluster. 
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Fig. 8. Gaps in temperature vs. age (upper panel) and 
metallicity (lower panel). See text for the metallicity un- 
certainties. 



